Shear induced breakup of droplets in a colloidal dispersion 
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We present numerical results for the breakup of a pair of colloidal particles enveloped by a droplet 
under shear flow. The smoothed profile method is used to accurately account for the hydrodynamic 
interactions between particles due to the host fluid. We observe that the critical capillary number, 
CaB, at which droplets breakup depends on a velocity ratio, E, defined as the ratio of the boundary 
shift velocity (that restores the droplet shape to a sphere) to the diffusive flux velocity in units 
of the particle radius a. For E < 10, Cas is independent of E, as is consistent with the regime 
studied by Taylor 0, [||. When E > 10, Can behaves as CaB = 2E~ X , which confirms Karam 
and Bellinger's hypothesis 0]. As a consequence, droplet breakup will occur when the time scale 
of droplet deformation 7 -1 is smaller than the diffusive time scale tr> = a 2 / Lr in units of a, where 
L is the diffusion constant and r is the 2nd order coefficient of the Ginzburg-Landau type free 
energy of the binary mixture. We emphasize that the breakup of droplet dispersed particles is not 
only governed by a balance of forces. We find that velocity competition is one of the important 
contributing factor. 

PACS numbers: 47.55.D-, 47.55.df, 82.20.Wt, 82.70.Dd 



A binary fluid system, such as a polymer blend, evolves 
under shear flow principally through the deformation and 
breakup of its constituent droplets Sheared binary 
fluids with dispersed colloidal particles are a topic of 
growing interest because they form the basis of many 
industrial products, for instance rubber products, plastic 
articles and some foods containing cornstarch or flour. 

Breakup of Newtonian droplets suspended in a New- 
tonian fluid was studied and first understood by Taylor 
U Q. He found, theoretically, in these purely viscous 
Newtonian systems that the deformation and breakup 
of droplets, in the absence of inert ial forces, is governed 
by the capillary number Ca = f] c Rj/a and the viscosity 
ratio p = r]d/r] c , where r] c is the viscosity of the contin- 
uous phase, R is the droplet radius, 7 is the shear rate, 
a is the interfacial tension and rjd is the viscosity of the 
droplet. The capillary number Ca = rjcRj/a is the ra- 
tio of the viscous stress exerted on a droplet by shear 
flow to the interfacial tension that attempts to restore 
the droplet's shape to a sphere. When the viscous stress 
is larger than the interfacial tension, the droplet becomes 
unstable and breaks. The capillary number at which this 
breakup occurs is known as the critical capillary number, 
CaB- Taylor suggested that CaB is a function of p. 

Taylor's result has been confirmed experimentally 0, 
and numerically 0, [Toj| . In the experimental works 
S SnSl, the behavior of Cas is in reasonable quantita- 
tive agreement. According to these results, Cas exhibits 
a minimum in the range 0.1 < p < 1. For 1 < p, within 
the highly viscous regime, Ca# increases with p until 
p « 4. For p > 4, droplet cannot be broken by shear. 
For 0.01 < p < 0.1, where the viscosity 77a is low, CaB 
monotonically increases with decreasing p. In the highly 



viscous regime, the rigidity of droplets increases with in- 
creasing p and finally the droplets behave almost like a 
rigid body. Thus, droplets cannot be broken by shear 
flow. In the low viscosity regime, Karam and Bellinger 
[3| hypothesized that Ca^'s behavior derives from the 
shape of deformed droplets prior to breakup. They ar- 
gued that a deformed droplet could be stabilized by the 
internally circulating fluid within it. Without internal 
circulation the deformed droplet is unstable and is likely 
to pinch off and divide due to interfacial tension. Inter- 
nal circulation counteracts pinching-off by building the 
pressure against dimples, pushing them outwards. This 
behavior is typical of softmatter systems in which slow 
and fast dynamics coexist. 

Evolution of the phase separation of binary fluid dis- 
persed colloidal particles has been investigated by many 
researchers [Ill4l7j. Introduction of colloidal particles 
strongly affects the phase structure. Neutral colloidal 



particles prefer to sit at the domain interface [ll|, \12 \. 



When the volume fraction of colloidal particles is suffi- 
ciently large, domain growth stops when the interfacial 
region is filled with particles. As a result, the final phase 
structure is a gel [12| . Colloidal particles introduced into 
a binary fluid suppress the hydrodynamic flow induced by 
interface motion if particles favor one of the fluid phases. 
As a result, phase separation slows down [13l-ll6j. With 
high concentrations of colloidal particles, attractive in- 
teractions, called wetting-induced depletion interactions, 
act between particles and thus the final structure be- 
comes a network of particles connected by one phase [TtJ . 

Droplets of a colloidal dispersion under shear have 
been investigated experimentally HI, [l9j and numerically 
20|, [H]. Inertia of the particles adsorbed to the droplet 
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interface causes the droplets to deform more strongly. As 
a result, droplets breakup more easily and Ca^ decreases 
in comparison to a Newtonian droplet [lo|. In two di- 
mensional simulations with low particle concentrations, 
the presence of dispersed particles or enveloped polymers 
serves to reduce the deformation of the droplet at in- 
termediate Ca. However, long polymers induce droplet 
breakup at high Ca [2l|. The mechanisms governing the 
breakup of droplets in the presence of dispersed particles 
are complex and have not been revealed completely until 
now. 

In a droplet of a binary fluid dispersed colloid under 
shear, internal circulation is strongly affected by the im- 
permeable surface of the colloidal particles. We hypothe- 
size that the presence of particles suppresses the diffusive 
flux within the droplet. This effect leads to droplet insta- 
bility. We analyze this effect and discover a new mecha- 
nism that governs droplet breakup. To accomplish this, 
we have devised a new methodology to directly simulate 
colloidal dispersions in a sheared binary fluid based on 
the smoothed profile method [lH, [23| . 

We extend the smoothed profile method 
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the dynamics of colloidal dispersions in Newtonian flu- 
ids to sheared binary fluids. In this method, boundaries 
between solid particles and solvents are replaced with 
a continuous interface by assuming a smoothed profile. 
This enables us to calculate the hydrodynamic interac- 
tions both efficiently and accurately, without neglecting 
many-body interactions. 

We employ the following free energy functional for a 
binary mixture containing particles: 



X (AV) 2 </>}, (1) 



where ip is a concentration field and (p is a smoothed pro- 
file function. The smoothed profile function < (p(r, t) < 
1 distinguishes between the fluid and particle domains, 
with <p = 1 in the particle domain and <p = in the fluid 
domain. These domains are separated by thin interfa- 
cial regions with a thickness characterized by £. The 
first term on the right-hand side of equation (pp) corre- 
sponds to a Ginzburg-Landau type mixing free energy 
for a binary fluid: f(ip) = -rip 2 /2 + uip 4 /4 + K\ Vip\ 2 /2, 
where r, \i and K are constants. The second term repre- 
sents the wetting interaction between the binary fluid and 
the particle surface, which is represented by |V^| 2 . W 
expresses the strength of the wetting interaction. With 
W > 0(W < 0), the phase with ip > 0(i/j < 0) favours the 
particle surface. The third term is introduced so that the 
concentration field inside each particle favors a surface 
phase, x 1S the coupling constant and Aip = ip — W/\W\ 
(W^O), ij) (W = 0). 

The time dependence of the concentration field ip and 
the fluid velocity u are given by 



chp_ 
dt 



-{u • V)tp + LV • {(I - nn) • V/x} 



(2) 



Pi 



du . __ 
- + (n.V)n 



-IpVjJL - V P + T]V 2 U + pi 0/ p + / S hear(3) 



under the incompressibility condition V • u = 0, where 
ji is 5F /5ip, I is the unit tensor, n is a unit vector field 
defined by — V0/|V0| and pf, 77, L and p(r,t) are the 
density, shear viscosity, diffusion constant and pressure 
field of the solvent, respectively. The operator (I — nn) 
directly assigns the no-penetration condition at interfa- 
cial regions. Note that, in this paper, we assume that the 
two phases have the same viscosity and thus the viscosity 
ratio p = 1. A body force <pf p is introduced to ensure 
rigidity of the particles and to enforce non-slip boundary 
conditions at the fluid/particle interface. Mathematical 
expressions for (j) and <pf p are detailed in previous pa- 
pers [iH, [HI . The external force /shear is introduced to 
maintain a linear shear with shear rate 7. This force is 
applied using an oblique coordinate transformation based 
on tensor analysis [24], . 

In the present study, we treat colloidal particles as 
beads of radius a. The motion of the ith bead is gov- 
erned by Newton's and Euler's equations of motion with: 
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where i^, Vi and are the position, translational veloc- 
ity, and rotational velocity of the beads, respectively. Mi 
and Ii are the mass and moment of inertia, and F^ 1 and 
iV 2 H are the hydrodynamic force and t orq ue exerted by 
the solvent on the beads, respectively [22|, [23|. Ff rep- 
resents the potential force due to direct inter-bead inter- 
actions, such as Coulombic or Lennard- Jones potentials. 
The truncated Lennard- Jones interaction is expressed in 
terms of U^j: 



ULj(r i: j) = < 



where r^- = 




(rij > 2% a), 



Ri — Rj\. The parameter e characterizes 



the strength of the interactions. 

Numerical simulations have been performed in three 
spatial dimensions with periodic boundary conditions. 
The lattice spacing A is set to unit length. The phase 
interfacial length I = yjKjr = 1 or 2 and the inter- 
facial tension a = \J~tK = rl. The unit of time is 
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FIG. 1: Schematic illustration of the rotation and breakup of 
a droplet. T is the tumbling period. The legend relates the 
color to the concentration, (a) Ca is sufficiently small and 
the droplet can rotate, (b) Ca is sufficiently large and the 
droplet breaks up before it can tumble. 



10' 

= (a/T))/(Lx/a) 
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given by pfA 2 /^, where 77 = 1 and pi = 1. The sys- 
tem size is L x x L y x L z = 64 x 64 x 32. The range 
of 7 is 2.0 x 10" 3 < 7 < 4.0 x 10" 3 and that of L is 
8.0 x 10 -2 < L < 1.8. The remaining parameters are as 
follows: r = u = K/l 2 = 0.01 - 1.0, x = 10W = IOO/777, 
a = 4 or 5, f = 2, e = 1, M» = 47ra 3 /3 and ft 0.064, 
where h is the time increment of a single simulation step. 

The initial state of the binary fluid is set according 
to the following procedure. At first, two particles are 
introduced at (32 — a, 32, 16) and (32 + a, 32, 16) into 
a homogenerous concentration of average concentration 
$ = -0.93 and r = u = K = 0.1, W = 0.04, X = 0.4, 
L = 1.6. We then iterate until steady state. Afterwards, 
we increment by a single time step. 

In the following simulations, the time dependence of 
the concentration field ip and the fluid velocity u are 
discretized using a de-aliased Fourier spectral scheme in 
space and an Euler scheme in time. To follow bead mo- 
tions, the position, velocity and angular velocity of the 
beads are integrated by the Adams-Bashforth scheme. 

Fig. [I] is a schematic illustration of the rotation and 
breakup of a droplet. When the interfacial tension is suf- 
ficiently large, as in (a), the droplet undergoes a tumbling 
motion similar to a chain [26|, although this tumbling 
has a limited lifetime. On the other hand, in (b), the 
interfacial tension is small and we expect the droplet to 
break-up before undergoing tumbling. 

We define the breakup shear rate 7b as the shear rate 
at which a droplet of a colloidal dispersion becomes un- 
able to tumble for a given a. The break up capillary 
number Ca# is defined as 



assuming that the droplet radius is roughly estimated by 
2a. 

To quantify the barrier to internal circulation due to 
the particles, we introduce the velocity ratio E which is a 
dimensionless number given by the ratio of the boundary 
shift velocity (that restores the droplet shape to a sphere) 
to the diffusive flux velocity in units of the particle radius 
a. Concretely speaking, E is expressed as 
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We plot the behavior of CaB as a function of E in Fig. 
[21 Were droplet breakup driven purely by the mecha- 
nism described in earlier works [ij-y, [5|-l9| , CaB would be 
independent of E. For E < 10, CaB is constant, namely 
CaB ~ 0.24. Previous works [lH3|, |5|-l9| reported a critical 
capillary number of CaB ~ 0.5 in the absence of particles 
and with equal droplet and solvent viscosities. Our value 
of CaB is smaller than the previously reported result due 
to the presence of colloidal particles near the interface. 

We confirm that the mechanism governing the breakup 
of droplet dispersed colloids in this regime is the same 
mechanism as described in Taylor's works [H, Q. Accord- 
ing to Taylor's work [l|, a critical capillary number CaB 
is estimated. We neglect the deviation of the droplet 
shape from a sphere. When the droplet radius is roughly 
approximated by 2a, the maximum pressure difference 
across the interface is Sp ~ 4t?7 and the Laplace pressure 
due to interfacial tension is Spi, = a/a. The Stokes force 
exerted on the colloidal particles is estimated as Gnrjadj 
where d is the distance between particles, approximated 
as d = 2a. Thus, the Stokes force per unit area, oTs, 
is approximately 67^07 /27ra 2 = 6777. Assuming that 
the droplet breaks when Sp + SF$ > 5pi, , we can predict 
CaB ~ 0.2. This estimate is in good agreement with our 
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FIG. 3: Cclb as a function of a 2A f/Lr. Closed circles (a = 
4,7 = 2.0 x 10~ 3 ,/ = 1). Open circles (a = 4,7 = 4.0 x 
10- 3 ,/ = 1). Closed triangles (a = 4,7 = 1.0 x 10" 3 ,/ = 2). 
Open triangles (a = 5,7 = 2.0 x 10~ 3 ,/ = 1). The dashed 
line corresponds to Cclb = 0.24. The solid line corresponds 
to a 2 7/Lr = 1.14. 



result Ca B ~ 0.24. 

However, when E 1 > 10, Ca B comes to depend on 
E. Although the interfacial tension exceeds the viscous 
stress, the droplet is broken by the shear flow in this 
regime. We reveal for the first time that Ca B can be 
described by 

Ca B = 1.80£T - 93 . (9) 

This equation implies that Ca B decreases with decreasing 
inner circulation due to the diffusion flux. This behavior 
confirms Karam and Bellinger's hypothesis [3|- 

Furthermore, Eq. (0 not only confirms Karam and 
Bellinger's hypothesis [3[ but also demonstrates a new 
source of instability originating from the competition 
between shear flow and diffusive flux. For droplet dis- 
persed colloidal particles, internal circulation, which cor- 
responds to a phase diffusive flux, is obstructed by the 
particle surface. The phase diffusive flux takes some 
amount of time, t& = a 2 /Lr, to overcome this obstruc- 
tion. When is larger than the time scale of droplet 
deformation 7 -1 , droplets can break up even when inter- 
facial tension exceeds the viscous stress. Therefore, we 
predict that the droplet is unstable when a 2 j/Lr > 1. 
On the basis of this, an instability condition can be de- 
rived as: 

Ca>2E~\ (10) 

which is in good qualitative agreement with Eq. ((9]). 

To verify the instability condition, Ca B is plotted as a 
function of a 2 j/Lr in Fig. 03 Fig. [3] shows that the 
droplet is stable for both Ca < 0.24 and a 2 j/Lr < 
1.14. Although the interfacial tension exceeds the vis- 
cous stress, the droplet is unstable and breaks up when 
a 2 j/Lr > 1.14. This result is in good agreement with 
our prediction, thus validating our instability condition. 



The breakup of pairs of colloidal particles enveloped 
by a droplet under shear flow was numerically calculated 
using the smoothed profile method extended to sheared 
binary fluids for 4.0 x 10 -2 < C a < 4.0 x 10 _1 and 2.0 x 
10" 1 < a 2 j/Lr < 2.0. Breakup occurs for 0.24 < Ca 
and 1.14 < a 2 ^/Lr. 

Droplet breakup at Ca ~ 0.24 occurs because the vis- 
cous stress exceeds the interfacial tension. 

At a 2 j/Lr w 1.14 breakup occurs because the time 
scale of droplet deformation 7 -1 is smaller than t& = 
a 2 / Lt. We emphasize that this instability is not derived 
from force competition. It is driven by velocity competi- 
tion. 
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